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We discuss elastic ed scattering beyond the Born approximation. The reaction amplitude contains 
six generalized form factors, but only three linearly independent combinations of them (generalized 
charge, quadrupole and magnetic form factors) contribute to the reaction cross section in second 
order perturbation theory. We examine the two-photon exchange and find that it includes two 
types of diagrams, where two virtual photons are interacting with the same nucleon and where the 
photons are interacting with different nucleons. It is shown that the two-photon exchange amplitude 
' is strongly connected with the deuteron wave function at short distances. 
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I. INTRODUCTION 



The study of electron scattering on the nucleon and the light nuclei provides a convenient tool to study the structure 
of strongly interacting systems. Due to the smallness of the fine structure constant a « — ^, one may expect that the 
Born approximation (one-photon exchange, OPE) should describe such processes with an accuracy of a few percent. 
Nevertheless, JLab polarization measurements of G P E (Q 2 ) / G P M (Q 2 ) [IHIJ], together with their theoretical analysis jH-Q], 
' show that higher order perturbative effects, such as two-photon exchange (TPE), can strongly affect some observables 
| of the elastic electron- nucleon scattering. 

Also for more complicated hadronic systems, like the deuteron, 3 He, 4 He, etc., TPE should contribute. Thus for 
precise studies of these nuclei a quantitative theoretical investigation of TPE effects is important; until now only a 
few estimates have been done of the contribution of TPE in. 

The aim of this paper is to estimate the TPE amplitude for ed-scattering in the framework of semi-relativistic 
calculations, with deuteron wave functions from "realistic" NN potentials. 

The paper is organized as follows. In Sect. |TT] we study the general structure of the reaction amplitude beyond 
^) OPE and define six independent generalized form factors which determine the amplitude. We show that only three 
linearly independent combinations of these generalized form factors contribute to the cross section in second order 
perturbation theory. We call the corresponding combinations of the form factors generalized charge, quadrupole and 
magnetic form factors. These generalized form factors are computed in Sect. [TTT1 Sect. ITVl contains numerical results 
?— i ' and a brief discussion. 



II. KINEMATICS AND DEFINITIONS 



The electron and deuteron momenta in the initial and final states of elastic ed scattering are denoted by k, k' 
and d, d! , respectively; q = k — k' is the transferred momentum; M and m are deuteron and nucleon masses; actual 
calculations will be done with m . 

All calculations arc done in the Breit frame, where the deuteron has the same energy Ed in the initial and final 
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state and moves along the z-direction (Fig. [T]). We get 



do = d' = E d = v/M 2 + Q 2 /4, 
d± = d' ± = 0, d 3 = -d' 3 = -Q/2, 

go = qi = 92 = 0, q 3 = Q, (1) 
k a = k' Q = E e , k ± = k' ± , k 3 = -k' 3 = Q/2, 

where Q is the modulus of the transferred momentum. For dcfinitcness we will assume that the transverse momentum 
of the electron is directed along the a;-axis 

ki = E e cos §, fc 2 = 0, k 3 = \Q = E e sin f . (2) 

In this frame the commonly used polarization parameter e can be expressed in terms of the electron scattering angle 

cos 2 | 

(3) 



1 + sin 2 § ' 



note that tg 2 § = (1 + ^tg 2 ^, where rj = 

The polarization vectors for the incoming and outgoing deuterons with spin z-projection A are denoted by £(\){d) 
and eix\(d'), respectively, 



c 



(d) = e (±1) (d') = -J±(0,±l,i,0) 



(±i 

e (0)( d ) = Yl l - ?' ' '^) . (4) 

mid 1 ) = ^(f.°> >^ 

The required electron electromagnetic current j M = M? l (fc')7 A 'u/ l (fc) is written as: 

jo = 2E e cos f , ji = -2E e , ^ 

j 2 = -2ihE e sin£, j 3 = 0. 

Because the electrons are ultra-relativistic, the helicities of the incoming and outgoing electrons are the same; its sign 
will be specified by h. 

Instead of the usual reaction amplitude M. it is useful to introduce the reduced amplitude Tyx,h by 

167TCt 

M = -yr^E e E d TyxM- (6) 

It follows from P and T invariance that, in the Breit frame, this amplitude must have the following properties 

Tyx-h = (-l) x ~ x 'T_y-X-h, from P invariance, ^ 
Tyx-h = T_ X -\';h, from T invariance. 
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The reaction amplitude is determined by six independent invariant amplitudes (form factors), which are specified by 
the following parametrization: 
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where lines and columns correspond to the allowing order: (A', A) = +1, 0, — 1, and 

1 i o = /i + /isin|/ 2) 0£_i=/3 + Asinf/4. (9) 

The form factors 0n, Qoo, fi, fi are complex functions of the two independent kinematical variables, for example 
Q 2 and 0. 

The relation of the amplitude Tx>x-,h to the invariant amplitudes G\, Gq used by other authors [3, [H| is given 
in Appendix [A] 

Later on the amplitude ([SJ will be expanded in a, and only terms of order zero and one will be kept. As shown in 
Eq. (|A8p . at zero order (OPE approximation) the amplitude written in terms of the charge, magnetic and quadruple 
form factors (Gc, Gm and Gq) becomes 



T (o) _ 



(Gc ~ frjGg) cos § -^/f G M (1 + hsin § ) 
§G M (l-frsin§) (G c + |^Gq)cos" 







^G M (l-/isinf) (G c - f »?Gq) co! 



VfGA/(l + /isin.! 



(10) 
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where the form factors Gc(Q 2 ), Gq(Q 2 ) and Gm(Q 2 ) are real and depend upon Q 2 only. 

Next the generalized electric, quadrupole and magnetic from factors Gc(Q 2 , 9), Gq(Q 2 ,&), Gm(Q 2 , 9), which repro- 
duce the spin structure of Eq. (|10[) . and the additional form factors gi(Q 2 ,8), g2(Q 2 ,9), 93(Q 2 ,9) are introduced as 
follows: 

Gn =Gc — 3VGQ, Goo = Gc + ^vGq, 
fi = Gm +gisin 2 f, fz = Gm - 9i, 

h =92, fi = g3- 

In OPE+TPE approximation the form factors can be written as: 

Gc = G c + SGc, Gq = G Q + 5Gq, 
Gm = G M + SGm, 

where d stands for the terms of order a; likewise, the form factors 51,2,3 are also proportional to a. 
By standard calculation one derives the differential cross section 

da 

dh 



cm 



in 
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(11) 
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where <7m is the Mott cross section and 

X,X',h 



X;h\ 
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[A(Q 2 ,9) + tg 2 ^B(Q 2 , 9)] + 0( 



= (l+sin 2 §) 



t\G E (Q 2 ,i 



:v\Gm(Q z , 
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with 



A(Q 2 



+ 0(a 2 ), 



\Gc(Q 2 ,t 



Iv 2 \g q (q 2 , 



Iv\Gm{Q 2 M\ 



B(Q 2 ,e) = l(i + v ) v \G M (Q 2 ,e)\ 2 , 



(15) 



G 2 E = \Gc{Q 2 M 2 + b 2 \QQ(Q 2 ^)\ 2 - 
The advantage of using the form factors Gc > Gq and Gm is that the expression for the cross section has the same form 
as the Roscnbluth formula; nevertheless the Rosenbluth separation of the structure functions A(Q 2 ,6) and B(Q 2 ,9) 
can no longer be done because they depend on two variables. 






FIG. 2: Two-photon exchange diagrams. The top diagrams correspond to the amplitudes M l v and M„, the bottom diagrams 
to the amplitudes (left) and (right). 

III. CALCULATION OF THE TWO-PHOTON EXCHANGE 

In what follows the contribution of meson exchange currents to the TPE amplitude will be neglected, and two types 
of TPE diagrams will be considered, where the virtual photons interact directly with the nucleons 



M 2 = M 1 +M U . 



(16) 



One of them, Ai 1 = A4 l p + M} n , corresponds to diagrams, where both photons interact with the same nucleon (Fig. [2l 
top). The other type, M u = M P l + corresponds to the diagrams, where the photons interact with different 

nucleons (Fig. [21 bottom). 

The deuteron structure will be described by the non-relativistic wave function 



= E CTllCT2 [yfi <52<^|1A) U (p)- (17) 

- E e ,M Y 2i (P) < 1 i^a 1 1M) (21£M 1 1 X)U 2 (p) 
x\N 1 a 1 ,N 2 a 2 ), 

where p is the internal momentum in the deuteron, |iVi<7i, N2O2) is the spin-isospin wave function of the two nucleons 
and (...I...) are Clebsh-Gordan coefficients. 



A. M l N diagram 



The TPE amplitude for a nucleon N has the following structure Q 



47TCV / 
M 2lN = -JyfU^^Uh (p' N (j' 



PN<J , 



(18) 



where H 1 ^ is the "effective hadron current" 



4m 



/??- 



(19) 
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In Eqs. (fT8|) and (fT9|) and pAr are the nucleon momenta, a and a' are the nucleon spin projections, |pjv f) and 
|pAr c') are the nucleon spinors, If = (k + k')/2, P = (p^ + pk)/2', AFf and AF^ 1 may be called corrections to the 
Dirac and Pauli form factors and F^ 1 is a new form factor. All the quantities AF-^ , AF^ 1 and F-^ are of order a. 
They are complex functions of two kinematical variables, e.g. Q 2 and v = APK. 
Because we employ a non-relativistic deuteron wave function we have to put 



(to, \d + p), 



Pp 
j 



Pn 



(to, \d-p), 
«(m, fd'-j?'), 



(20) 



where p and p' are the internal momenta in the deuteron. 
From (|18p it follows that the M l amplitude is given by 



(21) 



where X>^(A'A) is an effective deuteron current. The latter is derived in the same way as the deuteron current J' 1 in 
the impulse approximation 16!] with the nucleon current substituted by the effective hadron current (JH 



2?m(V, A) = — / dW(?+ y )( i/M + A). 

TO J 4 F 



(22) 



Here ^(p, A) and ^(p + ^(f, A') are wave functions of the deuteron in the initial and final states 1 . 

Later on we will need a non-relativistic reduction of matrix elements of the effective hadron current. Retaining the 
terms linear in the nucleon momentum one gets (see Appendix |B|) 



(pnct' 



H 



N 



p N a 



nxl ( 60S ~ cos § F 3 N ) Xa ee x lu%Xo 



2m 2 



Pnct 



H 



N 



PN <J 



(23) 



i (a x q) 6Q» 



where a are Pauli matrices and x&'-i Xo are Pauli spinors. 

The generalized nucleon electric and magnetic form factors are defined by (see Ref. 0]) 



dg% = AF? - tAF 



6G" = AFf + AFi 



N 



4m 



4to 2 



jv 

2 r 3 ' 



-F 



(24) 



Here r = « 4ry and v « mi? e . After substitution of Eqs. (|23|) and ([24]) into Eq. ([22]) the matrix elements of the 
effective deuteron current between the initial and final deuteron states become 

' 2E d (6g I c -h6g I g), if A = A' = ±1, 

2E d (sg l c + l v sg l Q ), if A = A' = 0, 

P°(A',A) = { -2i^x 
m 

Xy^cosKA'lJalA)^, if, A' - A = ±1, (25) 
, 0, if A' - A = ±2, 

V(\>, A) = i(\'\Jx r|A>§ fe, - ^ 



1 Here the additional multiplier (2m) 1 appears due to the additional multiplier 2m in the effective current [see Eq. l|23|l ] in comparison 
with the nucleon current (3) and (4) of Ref. [Trjl . 
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where J = (Ji, J 2 , J3) is an operator of the deutcron total angular momentum, and 

sg 1 c = 2sg% [i° 00 (Q) + i° 22 (Q)} , 



8Q l Q - 



n 

M 
111 



iio(Q) 



2V2 



iUQ) 



- } sg s E [i° 22 (Q) + iUQ)] 



2sg s 



M 

M 



W) - ¥UQ) + \ Pio(Q) + \ih(Q) 



(26) 



J- 3 = 2 Fn 



III 



Joo(Q)-i^2(Q)- 

In these expressions the notation l£l e (Q) = drjL (\Q r ) u i' (r)ue(r) was used, where jh{x) is a spherical Bessel 
function, ug(r) is the radial deuteron wave function for orbital momentum £ and Sg^ = ^(Sg^+Sg 1 ^), etc. Contracting 
the effective deuteron current with the electron current one arrives at 



-e — J3. 



(27) 



The corrections g 2 and g% are obviously vanishing for the first type of diagrams, g\ 
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B. M u diagrams 

In Ref. fl2| the contribution of the A4 U diagram was estimated within a non-relativistic approach with a Gaussian 
deuteron wave function. The present calculations are similar, but use a deuteron wave functions extracted from a 
"realistic" NN potentials; also, a modern paramctrization for the nucleon form factors has been adopted. 

The appropriate amplitude is given by the sum of two diagrams displayed at the bottom of Fig.[3J Ai n = Aip+Ai^^ 
where 

' M "* = / |}W ^'(A, A 2 ,Zl|^ . (28) 
Here p = \(j>^ — P^ 1 ) and p' = ^(p'^ 1 ' — p'^) are the relative momenta in the initial and final deuteron, 

~ P _ u h (k') (-iejf,,) itf, + n){-ie~f v )u h {k) 



t x =t F 
G(A lf A 2 ) = 



= (ie) 2 Tr{zd (A,) (p' (1) ,/ 2) ) +m)x (29) 
xr? (A 2 ) tfW +m)id,M(pM,pW)yW _ m ) x 
xf 2 ^(A 2 ) {fW-m)}, 

rr (A'A) ( A 2 jA 2 )=7r (A'A) (A 2 jA 2 ); 

D = [(p /(1) ) 2 -m 2 +i0][(p (1) ) 2 -to 2 +iO]x 
x [{p^f - m 2 + zO][(p' (2) ) 2 - to 2 + iO}. 

where P and X superscripts (subscripts) mean appropriate quantities related to diagrams with "parallel" photons (left 
bottom, Fig. 2) and "crossed' photons (right bottom, Fig. 2); 4 = ^ n p9j) we use the following notations, / is the 
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FIG. 3: Four types of poles taken into account in integration over dpo and dp' . 



4-momentum of the intermediate electron, Ai = k — I and A2 = I — k' are the 4- momenta of the virtual photons and 
k is an infinitesimal photon mass introduced in the photon propagators to regulate the infrared divergences; T^(Af) 
and T 2 (A|) are electromagnetic currents for the nucleon and anti-nucleon, in which the form factors are functions of 
A 2 and A 2 ,, respectively; d^\p^' ,p^ 2 ') and )(p / ( 1 ),p'( 2 )) are dpn vertex functions for the initial and final deuteron. 

In expression (|29|) for Tp^ X A ' moving along the nucleon loop (bold lines in the bottom diagrams in Fig. [2|: a line 
with an arrow opposite to the motion corresponds to a fcrmion propagator and a line with an arrow along the motion 
corresponds to an anti-fermion propagator. 

An infrared divergent term appears in A4 U when one photon is soft, Ai — > 0, A2 — > q or Ai — > q, A2 — > 0. It is 
canceled by radiative corrections which are not of interest in this paper. The configuration where each intermediate 
photon carries about half of the transferred momentum (hard-photon approximation) is emphasized, 



Ax 



A 2 ~^ 
2 



(30) 



In this case there is no infrared divergent terms. 

To relate the dpn vertex to the deuteron wave function with one of the nucleons on the mass shell one has to 
integrate over dpo and dp' . Four types of poles contribute to this integral (see Fig. [3J . What follows is a discussion 
of the contribution coming from the poles of the first diagram of Fig. [3] 



<£zr x 



d^pd^p' / P x 



{2nf 



Ti 



D 

pv(X A) 
P,X 



1 d 3 pd 3 p' 
~4 2E 1 2E' 2 (2ny 



(31) 



»o) (; 



(2Y 



iO 



where E± = J m 2 + (p— \q )" and E' 2 = J m 2 + {p' — jq )' 
One can use the expansion 



p'W+m = J2\P {1) <ri)(p {1) <ri\, 



(32) 



(in the last equation c means charge conjugated spinor) and define the wave functions of the initial and final deuteron 
by 



dW (»(!) D ( 2 )) 

^h^\p^)= J p f L , 

, s , s d( x ") (VW v'W) 
p'(i) 2 -m 2 + i0 

These wave functions are normalized by the condition 

/ 2^F Tr 0t (A) {p(1) ' p(2) ^ W {p(1) ' pi2) } = 1 

[and similarly for <fy^ ^ (p'^\ p'^ 2 ')], which comes from the requirement Gc(0) = 1. 



(33) 



(34) 
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Note that in general the nucleons N'^ and are not on-shell and at this step one cannot use expansions similar 
to for ^ + m and j) ( 2 ' — m. Nevertheless we will assume that the relative momenta in the initial and final 
deuteron are restricted by 

IpI ~ W'\ « Q. (35) 
This means that in all expansions one has to keep terms linear in p and p ' only and 



p^) = -\q±p, p'^) = \q±pi, (37) 



E\ y 2~\E d ±—- — , E' 12 K,\E d ± — — — , (36) 



Ai,2 = \q±8, S=^-^-(p> 3 + P3 ),p>-p). (38) 

One sees that in the framework of our approximation all nucleons become on-shell and one can use an expansion 
similar to (|3"2"j) for ft'* 1 ' + m and ft 2 * 1 — m. As a result all diagrams of Fig. |3]give the same contribution and 

„ e 2 d 3 pd 3 p' 
p ~ {2-k)*2E x 2E 2 

x <p' (2 M;c|^V' ( V)x (39) 

°" 1 i a 1 i a i i cr 2 

xip'^a^iA^p^a^ip^a^lp^^;^ x 
x(p^a 2 ;c\T^Al)\p'^a' 2 ;c). 

dT£ v is obtained by the exchange A 2 <-» A 2 in {39]). Using the fact that (jp Wai|</> A |p ^a 2 ; c) and (p' ( -^a 2 ;c\^ y \p' ( - 1 ^a' 1 ) 
are Lorentz invariants, one can substitute the non-relativistic deuteron wave functions (fl~7| instead of these wave 
functions: 

f d (40) 
L =(p'W* 2 -,c\^'\p'W*[) -+ (2n)VHl 



This substitution must be completed by the transformation of the current T 2 — > T 2 . In (j40|) 

(41) 



p'=(pi\P2',jf d Ps') = (p'x,jf d P 3 ' 

are the internal momenta in the deuteron rest frame. 

Expanding the current matrix elements (p'W er'|rf (Af)|p W a) = 2mx^rf (Af)xa in terms of Pauli spinors one 
gets 

„ _ 64a 2 (4n) 2 E d f d 3 pd 3 p' 

M Q e J 1^F X 

x [rff lAl (iQ 2 )r 2 „(lQ 2 ) + $A+ (42) 

+p"B + <D(p 2 7 p' 2 )] *(A,p), 
where A and £? are some vectors and 

The integrals J d 3 p p^>(X,p) and J" d 3 p' p'*S>(\' ,p' ) obviously vanish after angular integration and one arrives at 

MU .- ^y- ^y, (44) 
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where 9Jt^ A = ^*x,{0)r lti (jQ 2 )r2v{jQ 2 )ip\(0); ip\>(0), ip\(0) and il>\>(0) are the deuteron wave functions in coordinate 
space at r = and 

f° (±Q 2 ) = G% (IQ 2 ) , 

r,QQ 2 ) = ^x^(iQ 2 ) (45) 

(suffix k = 1,2 enumerates the nucleons). 
For further calculations it is useful to introduce "plus" and "minus" components of the tensors according to A± = 

\{A\ ± 1A2). The contraction of the lepton and deuteron tensors becomes 



where 



and 



r 00 = 8Ej cos f , 

r 0+ = -2V2E 2 (2 - sin 2 § - /i sin |) , 
to- = -2\/2X 2 (2 - sin 2 f + /isin |) , 
r ++ = 4£ 2 cos § (1 - h sin f ) , 

r = cos § (1 + /isinf) , 

r_ + = 4£f cos f 



mil = ^o^ 1 = - OT o- = -^o-° = 



2tt 

47T 



with the abbreviations 



C= K(r)] 2 



r=0 

G BB = G|(iQ 2 )Gl(iQ 2 ), 
Gmm = G p M (\Q 2 )G 1 l I (jQ 2 ), 
Gem = \ [GUlQ^GUlQ 2 ) + G P M (\Q 2 )G%(\Q 2 )] 
Finally we get the following amplitudes 

Til = * cos §G EE , 

T oo = KC0S I {Gee - vGmm) , 

T w,h = -x^Gem (2 - sin 2 § + hem §) , 

^l-i./i = xvGmm cos § (l + /i sin I) , 

where 

l28aCE e 



(46) 



(47) 



-^V^Gbm, (48) 



(49) 



(50) 
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and one arrives at 

8Qc = x (Gee - \vGmm) 



xr ii _ 2hGem 

* Gm l + sin 2 !' (51) 



n _ kGem cos 2 § 

91 ~ l + sin 2 § ' 

92 = 93 = cos | G M M- 

One should note that the approximation (j42|) is not valid at — > (or, equivalcntly, e — > 1). Indeed, the denominator 
of the electron propagator 

^O 2 + ^(p* +P'J - 2cosfi? e fe -p' x )Q* 

contains products of E e and components of the internal momenta. From ^ it follows that E e — > oo when 9 — > 
and the factor C should be changed to 

C^S = 

1 f d 3 pd 3 p'Uo(p)U (p'^ (52) 

The amplitudes P5j) with the substitution (|5^|) in (|50|) coincide with the results of Ref . [l2[ . 
To evaluate the integral above one can use the integral representation for the denominator 

i 

dre l{a+M)T (53) 



a + iO 

and reduce (|5^|) to a one-dimensional integral 



roc j 

5=-|Q 2 / ^% 2 (y), 



where ?/ = T-E e ^/ 4 cos 2 f + . Changing the variable in (|54j) one gets 



S = -if r%e^ul{y), 

Jo v 



(54) 



(55) 
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TABLE I: Uq(0) for some popular potentials. 



u '(0), fm 6,1 


Potential 


Ref. 


1.1978 ■ KT 1 


Paris 


[18] 


3.1035 • KT 1 


CD-Bonn 


ris] 


2.6860 • 10" 1 


Nijm I 


[20] 


2.6730 ■ 10 -2 


Nijm II 


[20] 


3.1571 • 10" 1 


Nijm 93 


[20] 


5.8334 • 10" 2 


Reid 93 


[20] 



where 



O 2 

f= , = (56) 

For the standard parametrization of the wave function 

u (y) = ^2c n e- a ^, with ]Tc„ = (57) 

n n 

we obtain ( Appendix |C|) 

S = -if ^2 2J c„c m (a„ + a m - i/) ln(e*n + «m - £/)• (58) 

n m 

In (|55[) the exponent reduces to 1 in the limit E e — > oo for fixed Q and 

9?eT# A ~ sin §, 3m^! A -> consi, (59) 

i.e., at the limit the TPE does not contribute to the cross-section in the next order of the a-expansion. 

The e and Q 2 dependence of S is displayed in Fig. 21 One sees that the S factor depends strongly on the NN 
potential and in any case it is very different from the constant value C = [u' Q (0)] 2 . The reason is as follows: from (|58|) 
one gets that S — > C at the formal limit 

(60) 

But from expression (f56|) it follows that / — > iMsin| when Q 2 —> oo and the condition ([60)) cannot be fulfilled at 
any Q 2 . Note that a similar situation takes place in the evaluation of the so-called triangle diagram in pd backward 



scattering [17 



We have also studied the deuteron D-wave contribution to the S factor and found that it contributes less than 10%. 



IV. NUMERICAL RESULTS AND DISCUSSION 



Figs. [S][5] display the e-behavior of the TPE corrections SQc/Ge, $Gq/Ge, 5Qm/Ge and Qx^IGe calculated with 
the deuteron wave function for the CD-Bonn and the Paris potentials. The form factor Ge(Q 2 ) was calculated in the 
framework of the impulse approximation. 

In the present calculations of the TPE correction in Ai 1 the amplitudes AF/2 and F3 from the theoretical calcu- 
lations of Ref. [6j are used. At Q 2 < 6 GeV 2 they are practically independent of the parametrization of the nucleon 
form factor. For A4 U we use the following parametrization of the nucleon form factor: 

• for the magnetic form factors of the proton and neutron — dipolc parametrization 

G P M (Q 2 ) = V P G D (Q 2 ), G n M (Q 2 ) = (i n G D (Q 2 ), (61) 



where G D (Q 2 ) = (l + Q 2 /0.7l) 
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FIG. 6: Same as Fig. H for Q 2 =2 GeV 2 . 
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FIG. 7: Same as Fig. for Q 2 =3 GeV 2 

0.1 




FIG. 8: Ratio of gi = gs to Ge- Dashed, dot-dashed and solid curves are for Q 2 =1, 2 and 3 GeV 2 , respectively. Left and 
right panels are for the CD-Bonn and Paris potentials, respectively. 



• the electric form factors of the proton and neutron were taken from the parametrization of the JLab data (see 

M) 

G P E (Q 2 ) = (1.0587 - 0.14265Q 2 ) G D {Q 2 ), (62) 
and and so-called Galster parametrization [22j , respectively. 

One sees that two-photon exchange may give a large contribution to the elastic ed scattering; although caution is 
required as these estimates have large uncertainties. The most important source of uncertainty comes from the S 
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factor which is determined by the short range part of the deuteron wave function. The last quantity is very poorly 
known (see, e.g., Table HJ). Of course, besides NN degrees of freedom, non-nucleon (quark) degrees of freedom should 
also be taken into account in this region and one may expect that in the framework of more realistic estimates the 
two-photon corrections may be smaller. The implication is that the experimental study of two-photon exchange in 
elastic ed scattering at Q 2 ~ few GeV 2 can give important information about the deuteron structure at short distances. 

In summary, we estimated the two-photon exchange amplitude in elastic ed scattering. There are six independent 
form factors which determine this amplitude, but only three of them contribute to the cross section in second order 
perturbation theory 

There are two types of two-photon exchange diagrams. For the first type two intermediate photons interact with 
the same nucleoli. For the second type the intermediate photons interact with different nucleons. 

We show that the two-photon exchange amplitude is strongly connected with the deuteron structure at short 
distances. 
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Appendix A: General expression for the ed —> ed amplitude in the Breit frame 

From the invariance under Lorentz transformations and space and time inversions it follows that the amplitude of 
the elastic scattering of a spin-^ particle (the electron) off a spin-1 particle (the deuteron) has 9 invariant amplitudes 
(form factors). Usually for ultra- relativistic electrons the mass can be neglected and the electron helicity is conserved. 
In this case the number of form factors is reduced to six and the amplitude has the general form [3, HH 

Ana ,„ Ana . 

<T Q 

where j M is the e.m. current for the electron and J M is an "effective current" for the deuteron 



m = — = -=3-r VAih> (Ai) 



-{G x (e'*-e) (d + d')^ + 

+ G 2 [(e'*-g) e,-e';{e-q)]- 

(e'* • g) (e • q) , 
~ 2M 2 { ^ 

(e'* ■ K) (e-K) , (A2) 

+ G4± ^ L {d + d%+ 

j* 



G 5 [(e'* • K) e, + ^ (e • K)] + 

(eT-K) (e-q)-(e'*-q) (e • K ) 
2M 2 



Here K = k + k' . The form factors Gi, Gg are complex functions of two variables, e.g., Q 2 and 9. 
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In the Brcit frame one simply finds 



£(±) ' Q = e (±) '9 = 0, 



£(o) ' Q = e (0 ) ' 9 



E d Q 
M ' 



£ (±) ' e (±) = -!> £( T )-e(±)=0, 
£(±) ' £(0) = £(±) ' e'(o) = °> 

12 



6 (o) ' e (o) = - I 



Q \ 

2M 2 J ' 

£(±) ' K = £(±) ' K = ±^ E e 



COS 



2 ■ 



-(0) 



K 



-£ ( o) ' K 



E e Q 
M 



and 



E e Q a 



£(0) ' J = 
£(±) ■ J = \/~2E e (±l — /isin|), 
£(±) ' J = V2^ e (±l + /isin§), 
+ • j = 4£ e £ d cos§. 



£(o) ' J = cos ' 



A/ 



From Eqs. (|A"Tj) - (|A"4|) it follows that 



2 



Gi - 



M 



., 6 Ep 
cos' -G 4 - -f G 5 



COS : 



Too,/* = [(1 + 277) d - 2 V G 2 + 2(1 + 77)7^3- 
M 



. E e E d 

T]&4 — 4 ; f0 Gg 



10, /l 



= x/2 



-T, 



01, -ft 



M 2 



-T_io,-/i 



Q 



-E d G 2 + 2^$- cos 2 |G 4 + 



M 2 

+ E e (l + cos|)G 5 + 2^ 
+/isin§(-£; d G 2 + £; e G 5 )] , 



cos % G 6 + 



■j^T C0S 2 G 4 



1 + A sin I G £ 



(A3) 



(A4) 



(A5) 



In the Born approximation the form factors G4 , G5 and Gg vanish and the form factors G\ , G2 and G3 become real 

,(0), 



functions Gf\ and G 3 °^ of one variable Q 2 , that is 



(A6) 



G : =Gl 0) (Q 2 ) + O(a), G 2 

G 3 = G 3 0) (Q 2 ) + 0(a), G 4 - G 5 - G 6 - a. 

Commonly the charge, magnetic and quadrupole form factors, Gc(Q 2 ), Gm(Q 2 ) and Gq(Q 2 ) are used instead of the 
form factors G^ , G^ and Gg '' . They are connected by 



G i 0) (Q 2 )=G c (Q 2 )-|7 7 G Q (Q 2 ), 
G^ 0) (Q 2 )=G M (Q 2 ), 



G 3 0) (Q 2 ) 



l + r? 



[-G c (Q 2 ) + Gm(0 2 )+ 



(A7) 
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and Eqs. (|A5|) arc reduced to 



T- 



(0) 
11, h 



T 



(0) 

■1—1, h 



[Gc(Q 2 ) - §??G Q (Q 2 )] cos ■ 



[Gc(Q 2 ) + |^G Q (Q 2 )] cos. 



T- 



(0) 
10, h 



-T, 



(o) 



r, 



(0) 



01, -h ~ ''O-l.i 



-T 



(o) 

10,-/i 



= -a/^Gm(Q 2 )(1 + ft sin I) 



(0) 



1-1, h 



T 



(0) 

-11, -ft, 



0. 



(A8) 



Appendix B: Non-relativistic reduction of the effective hadron current 

In the Breit frame K = E e (l, cos |, 0, 0) and 

(jT' ( *V \K\ p W a) « x t, { 2m Ai? 1N + 

+ 7T- (-Q + 2ie 3nm p n a m ) AF 2N + 
2m 

E 1 

+— [2m - cos §(2^ + iQa 2 )] F 3N \ X a, 

p'Wa'm\pW a) *2xl [(p a - 
-iQe a3n a n )AF 1N + iQe 3an a n AF. 
p'^a'\H 3 N \p^a)^ 



"2N 



where a = 1, 2. 

Note that up to terms of order O 



(Bl) 



(B2) 



(B3) 



50 the amplitudes A-Fuv, Ai^jv an d are independent on the 

nuclcon momenta. In this approximation the terms proportional to p±_ will vanish after integration in Eq. (|22[) and 
Eqs. (|B2]) and (|B3j) become 



2mAF 1JV - ^-AF 2N - 
2 m 



'2m - £Q cos fcr 2 ) F 3N 



= xl \ 2mAg E - i^- cos f a 2 F 3Af ) 

p'^a'\H%\p^a)^ 
« 2 4 ge 3a "(AFiAr + AF 2W ) X >" Xo 



(B4) 



(B5) 



Let us consider the integral 



Appendix C: 



J = 



dy 

o y 2 



,ifv 



where / is a real constant. With the expression (f5T|) for uo(y) the integral becomes a series 



dy 

o y 2 



,(if-ot„-a m )y 



(CI) 



(C2) 



with all terms divergent. We regularize the expressions defining 



/ = lim h , 



where 

dy 



„2-£ 

o y 



J* [u (y)} 



2 



= ^2^}2cnC m (a n + a m ~ if) 1 T(-l + e). 

n m 

Expanding the T-function near the pole 



r(-l + e) = -i +7-1 -0(e) 

e 



and taking into account the constraint J2 n c n = one gets 

I = lim y~] y~] c n c m (a n + a m - if) ln(a„ + a m - if) + 0(e) 

n m 

= ^2^2 c n c m (a n + a m - if) ln(o! n + a m - if). 
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